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Abstract 

Global  uniqueness  of  inverse  problem  for  a  hyperbolic  PDE  with  the  unknown 
potential  and  a  single  incident  plane  wave  is  a  well  known  long  standing  open  question. 

This  question  is  addressed  for  the  case  when  derivatives  with  respect  to  variables 
orthogonal  to  this  plane  wave  are  expressed  via  finite  differences. 

This  preprint  is  avalaible  on-line  at  http://www.ma.utexas.edu/mp_arc,  2011,  posting 
date  May  3,  2011. 

1  Introduction 

For  the  sake  of  brevity  we  consider  here  only  the  3-D  case.  All  theorems  below  have  almost 
identical  formulations  and  proofs  for  the  n— D  case  with  n  >  2.  Below  x  =  ( x,y,z )  G  M3. 
Let  the  function  a  G  C2  (M3)  and  is  bounded  in  M3  together  with  its  derivatives.  Consider 
the  Cauchy  problem 

utt  =  A u  +  a  (x)  u,  (x,  t)  G  I3  x  (0,  T) ,  (1.1) 

u  (x,  0)  =  0,  ut(x,  0)  =  5  (z) .  (1.2) 

Conditions  (1.1),  (1.2)  mean  that  the  wave  field  u  is  initialized  by  the  plane  wave  at  the 
plane  =  0}  .  This  plane  wave  propagates  along  the  z— axis.  Let  Ax,  Ay,  T  =  const.  >  0. 
Define  the  strip  G  as 

G  =  (x  :  x  G  (0,  Ax) ,  y  G  (0,  Ay)}  ,  Gt  —  G  x  (0,  T) , 

ST  =  {z  =  0,z  G  (0,  A^)  ,y  G  (0,  Ay)}  x  (0,T) . 

Inverse  Problem  1  (IP1).  Assume  that  the  function  a  (x)  is  unknown  in  G.  Determine 
the  coefficient  a  (x)  for  igG,  assuming  that  the  following  two  functions  r  (x,  t ) ,  s  (x,  t )  are 
given 

u  |St=  r  (x,  t) ,  uz  |St=  s  (x,  t) .  (1.3) 
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Uniqueness  for  IP1  is  a  well  known  long  standing  open  question.  IP1  is  a  Coefficient 
Inverse  Problem  (CIP)  with  single  measurement  data.  We  use  the  term  “single  measure¬ 
ment”  since  only  a  single  direction  of  the  incident  plane  wave  is  used.  Also,  (1.3)  is  the 
backscattering  data.  The  main  challenge  is  the  single  measurement,  not  the  backscattering. 
Assume  for  a  moment  that  the  Fourier  transform  with  respect  to  t  can  be  applied  to  the 
function  u,  and  the  resulting  function  p  (x,  k)  satisfies  the  equation 

Ap  +  k2p  +  a  (x)  p  =  —  5  (z)  (1.4) 

as  well  as  radiation  conditions  at  the  infinity.  Then  one  can  derive  from  (1.3)  and  (1.4)  the 
statement  of  a  well  known  CIP  in  the  frequency  domain.  In  accordance  with  the  paper  [7], 
the  question  about  the  uniqueness  of  this  CIP  was  posed  by  I.M.  Gelfand  in  1954.  In  [7] 
this  question  was  addressed  for  the  case  of  infinitely  many  measurements.  We  also  mention 
works  [3,4],  where  a  CIP  for  a  more  general  hyperbolic  equation  was  considered  and  solution 
was  constructed  for  the  case  of  infinitely  many  measurements. 

As  to  the  case  of  CIPs  with  single  measurement  data  in  n  —  D  ( n  >  2),  currently  only  one 
class  of  uniqueness  theorems  for  them  is  known.  All  these  theorems  are  proven  under  the 
assumption  that  at  least  one  initial  condition  does  not  vanish  in  the  entire  domain  of  interest 
fh  They  are  proven  via  Carleman  estimates.  The  idea  of  application  of  Carleman  estimates 
to  proofs  of  uniqueness  results  for  CIPs  with  single  measurement  data  was  originated  in  [9] 
with  many  follow  up  publications  of  many  authors,  see,  e.g.  [2,5,10-15]  and  references  cited 
there  for  some  of  these  works;  the  most  recent  survey  of  these  results  can  be  found  in  [20]. 

In  this  paper  the  uniqueness  question  is  addressed  for  a  closely  related  inverse  problem. 
Specifically  we  assume  that  derivatives  with  respect  to  (x,  y )  are  written  via  finite  differences 
with  the  grid  step  sizes  ( hi,h2 )  ■  Numbers  hi,h2  do  not  tend  to  zero.  However,  derivatives 
with  respect  to  z,  t  are  written  in  the  usual  form.  Under  these  conditions  we  prove  uniqueness. 

First,  we  prove  a  Carleman  estimate,  which  contains  a  new  feature  of  the  positivity  of  a 
certain  integral  over  the  characteristic  curve.  Next,  a  combination  of  this  new  feature  with 
a  modification  of  the  method  of  [9]  leads  to  the  desired  uniqueness  result.  It  is  worthy  to 
mention  here  that  discrete  Carleman  estimates  are  attracting  an  interest  nowadays,  see,  e.g. 
[8].  However,  they  were  not  yet  used  for  proofs  of  uniqueness  of  discrete  CIPs.  A  discrete 
Carleman  estimate  is  not  used  here. 

Both  classical  forward  problems  for  PDEs  as  well  as  CIPs  are  often  solved  numerically 
by  the  Finite  Difference  Method  (FDM).  However,  there  is  a  substantial  difference  between 
forward  and  inverse  problems.  Indeed,  classical  forward  problems  are  well-posed.  Because 
of  this,  an  important  topic  of  study  in  this  case  is  the  convergence  of  the  solution  obtained 
by  the  FDM  to  the  actual  solution  when  the  grid  step  size  tends  to  zero  [19].  However, 
the  computational  experience  of  the  author  [16,17]  shows  that,  because  of  the  ill-posedness, 
there  is  of  a  little  help  to  investigate  the  convergence  of  FDM-based  numerical  methods  for 
the  case  when  the  spatial  step  size  hsp  tends  to  zero.  Unlike  well-posed  problems,  in  the 
ill-posed  case  hsp  should  usually  be  limited  from  the  below.  The  reason  of  this  limitation 
is  that  hsp  serves  as  an  implicit  regularization  parameter  in  most  cases.  Furthermore,  it 
is  an  important  observation  of  numerical  studies  that  hsp  usually  cannot  be  significantly 


2 


decreased.  For  example,  while  accurate  computational  results  were  obtained  in  [16]  for  the 
case  hsp  =  0.01,  an  attempt  to  decrease  hsp  to  0.005  led  to  a  significant  degradation  of 
images,  see  Remark  7.1  in  [16].  A  similar  observation  took  place  in  subsection  5.2  of  [17]. 
Therefore,  hsp  should  be  bounded  from  the  below  by  a  positive  constant,  which  is  the  case 
of  the  current  work. 

As  to  the  application  of  IP1,  consider  the  following  analog  of  the  Laplace  transform  [18] 

u  (x,t)  =  2J-t 3/2  J  u(x,r)rexp  dr  :=  Lu.  (1.4) 

o 

It  is  well  known  that  [18] 


c(x)Ut  =  A  U  +  a(x.)U,  (1.5) 

U  (x,  0)  =  5  (. z ) . 

The  transform  (1.4)  is  one-to-one.  Thus,  if  the  uniqueness  of  IP1  would  be  proven,  then  the 
uniqueness  of  a  similar  inverse  problem  for  equation  (1.5)  would  be  proven  as  well,  and  vice 
versa.  On  the  other  hand,  equation  (1.5)  governs  light  propagation  in  a  turbid  medium,  such 
as,  e.g.  biological  medium.  In  this  case  U  (x,  t)  is  the  light  intensity  and  (—a  (x))  >  0  is  the 
absorption  coefficient  of  light  [1],  The  absorption  coefficient  contains  an  information  about 
the  blood  content,  which  is  important  for  medical  imaging. 


2  Formulations  Of  Results 


Consider  partitions  of  intervals  x  G  (0,  Ax) ,  y  e  (0,  Ay)  in  small  subintervals  with  step  sizes 
hi  and  h2  respectively, 


0  =  x0  <  Xi  <  ...  <  xNl  =  Ax,  0  =  y0  <  V!  <  ...  <  Vn2  =  An  A1) 

Xi  -  Xi_ i  =  hr,  yj  -  yj-i  =  h2,  h  :=  (hi,  h2) ,  h0  =  min  (hi,  h2)  ;Ni,N2>2.  (2.2) 

Hence,  we  have  obtained  the  grid  G}x  =  {(x,  y)  :  x  —  ihi,  y  =  i/i2}|^wo0)  •  Consider  a  vector 
function  fh  (z,  t )  defined  on  this  grid,  fh  (z,  t )  =  {fid  (z,  .  For  two  vector  functions 

fh  (z,  t ) ,  gh  (z,  t )  define 


gh  (z,t)  fh  (z,t) 


{ hj  ( z,t )} 


(JVi,JV2) 
(*d)=(o,o)  > 


ki,j  (z,  t ) 


9i,j  {z,  t)  ■  fi,j  (z,  t ) . 


Denote 


(. fhM)2 


(NUN2) 

fiAz A)- 


(iJ)=(  0,0) 


(2.3) 
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Let  B  be  a  Banach  space  of  functions  depending  on  (z,  t ) .  Denote  Bh  the  space  of  above 
vector  functions  fh  (z,t)  such  that  ffj  (z,  t)  e  B,\/(i,j)  G  [0,  Ni]  x  [0, -/V2] .  We  define  the 
norm  |  \fh  (z,  i)||B  in  this  space  as 


(JVi,JV2) 

E  ii/«ii2b 

diJM  0,0) 


1/2 


We  dehne  hnite  difference  second  derivatives  dl  hfh  ( z,t )  and  dy  hfh  ( z,t )  with  respect  to 
x  and  y  respectively  in  the  usual  way  as 


dl,Jh  ( z,t )  =  (^^)}(^)=(o,0)  ’  dy^k  ^  =  (z^)}pj)= 


{NuN2) 


(0,0)  ’ 


f)2  f ■  ■ 


GO 


{z,  t )  -  2 fij  (z,  t )  +  /i+ij  (2,  t)  if  i  ±  0,  X  i-  Nu 
fij  (z,  t )  -  2/i+1j  (z,  t)  +  fi+2,j  (z,  t)  if  i  =  0, 
fi,j  (z,t)~  Vi-lj  {z,  t)  +  fi—2,j  ( Z,t )  if  i  =  w 


and  similarly  for  d2hfh  ( z,t )  .  Hence,  if  a  function  g  ( x,y,z,t )  has  continuous  derivatives  up 
to  the  fourth  order  with  respect  to  x,  then  d2  hgij  (z,  t )  approximates  gxx  (x,  y,  z,  t )  at  the 
point  (x,  y)  =  (ihi,jh2)  with  the  accuracy  O  (hf) ,  hi  — >  0  in  the  case  when  ih\  ^  0,  Ax.  And 
it  approximates  with  the  accuracy  O  (hi) ,  hi  — >  0  in  the  case  when  ih\  =  0,  Ax.  Similarly 
for  the  //-derivative.  Next,  we  dehne  the  hnite  difference  Laplace  operator  as 


A hfi,j  (Z,t)  ■  (A  A  A h,x,yfi,j  ( Z,t )  , 

A h,x,yfi,j  {z,  t)  .  d Xthfi,j  (z,  A  9ythfi,j  (z,  t)  j 


Ahfh  ( z,t ) 


=  {Aft/ij  (z,t)} 


(iVi,JV2) 

(m)=(o,o) 


=  ( d2Jij{z,t )) 


(NUN2) 

(iJ)=(  0,0) 


A  {A h,x,yfi,j  (z,  t'j } 


(JVi,JV2) 

(*J)=(0,0) 


=  d2Jh(z,t)  +  Ah,xJh(z,t). 


Dehne  ah  (z)  :=  { a (^)}[fjL)'=(o  o)  ■  Rewrite  the  problem  (1.1),  (1.2)  in  the  hnite  difference 
form  as 

<  =  A  huh  +  ah  (z)  uh ,  (z,  t)  eRx  (0,  T) ,  (2.4) 

Mh(z,0)  =  0,4(z,0)  =  <J(z)T  (2.5) 

where  the  product  cA  (z)  uh  is  understood  as  in  (2.3). 

Inverse  Problem  2  (IP2).  Let  the  vector  function  uh  (z,t)  be  the  solution  of  the 
problem  (2.4),  (2.5).  Determine  the  vector  function  ah  (z)  assuming  that  the  following  two 
vector  functions  rh  (t) ,  sh  (t)  are  given 


uk(0,  t)  =  rh(t),  uj(0 ,t)  =  sh(t),  t  e  (0,  T)  . 


(2.6) 
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Theorems  1  and  2  are  main  results  of  this  paper. 

Theorem  1.  Let  the  vector  function  ah  (z)  G  C1  (M)  and  is  bounded  in  M.  Then  there 
exists  unique  solution  of  the  forward  problem  (2.4),  (2.5)  of  the  form 

Ui,j(z,t)  =  ^H(t-  \z\)  +uid(z,t)  ,(i,j)  G  [0,  Ni]  x  [0,  N2\ ,  (2.7) 

where  H  (z)  is  the  Heaviside  function  and  the  function  ul:J  is  such  that 

uid  G  C3  (t  >  |z|) ,  Uij  ( z ,  t)  =  0  for  t  G  (0,  |z|] .  (2.8) 

Theorem  2.  Suppose  that  there  exists  two  pairs  of  vector  functions  [u\  (z,  t ) ,  a J  (z))  , 
(u%(z,t)  (z)}  such  that  G  C1  (M)  and  vector  functions  are  solutions  of 

the  problem  (2.4),  (2.5)  with  a \  and  a %  respectively  of  the  form  (2.7),  (2.8).  In  addition, 
assume  that  both  vector  functions  Ui,u %  satisfy  the  same  conditions  (2.6).  Let  R  >  0  be  an 
arbitrary  number  and  T  >  2 R.  Then  af  (z)  =  a %  (z)  for  \z\  <  R  and  u ±  ( z,t )  =  u %  ( z,t )  for 

(z,  t)G  {(l-|)  \z\  +  *t<R,t>0 

3  Proofs  of  Theorems  1  and  2 


3.1  Proof  of  Theorem  1 


Denoting  temporarily  /jj  (z,t)  =  Ahxyuh]  +  rewrite  (2.4),  (2.5)  as 

d"  fi,j  (A  ) 

Uij  (z,  0)  =  0,  dtUij  (z,  0)  =  6  (z) . 

Hence,  D’Alembert  formula  implies  that  for  (i,j)  G  [0,  Ni]  x  [0,  AT2] 

t  t—T+Z 


Uij(z,t)  =  ]-H(t-  \z\)  +  \ 


(It  I  (Ah  X  yUi  j  +  di  jiii  j )  (£,  t)  d^. 


0  T  —  t~\~Z 


(3.1) 

(3.2) 


(3.3) 


In  (3.3)  the  integration  is  carried  out  over  the  triangle  A  (z,t)  in  the  (£,  t)  —  plane,  where 
A  (z,t)  has  vertices  at  (^i,Ti)  =  (z  —  t,  0) ,  (£2)  t2)  =  ( z,  t )  and  (£3,t3)  =  (z  +  t ,  0) .  The 
set  of  equations  (3.3)  considered  for  (i,j)  G  [0,  Nf  x  [0,  N2\  is  a  linear  Volterra  system  of 
coupled  integral  equations.  It  can  be  solved  iteratively  as 


u 


(n) 


(z,t)  = 


\H{t 


kl)  +  i 


t  t—T+Z 

fdT  / 

0  T  —  t+Z 


A(n—  1)  . 

h,x,yuij  +  a 


ijuinj  1})  (f  > T )  d(- 
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Let  rnaxjj  supn  \ahJ  (z)  <  M,M  =  const.  >  0.  The  standard  technique  for  Volterra  equa¬ 
tions  leads  to  the  following  estimate 


u 


(n) 

ij 


[z,t) 


sE 

n= 0 


0 cty 

n\ 


,2  g  R,t  >  o 


(3.4) 


where  the  constant  C  =  C  ( h ,  M)  .  Hence,  there  exists  a  solution  of  the  integral  equation 
(3.3)  such  that  this  solution  is  bounded  for  any  t.  Let  utJ  ( z,t )  =  utj  ( z,t )  —  H  (t  —  \z\)  /2. 
Then  (3.3)  implies  that  the  function  Ttij  ( z,t )  G  C1  (t  >  |z|) .  Differentiating  equation  for 
uid  ( z,t ) ,  one  obtains  that  ut^  ( z,t )  G  C3  (t  >  |^|) ,  which  is  the  same  as  (2.8).  Thus,  the  set 
of  functions  u^j  ( z,t ) ,  which  solve  (3.3),  satisfy  conditions  (3.1),  (3.2),  (2.7)  and  (2.8). 

Next,  let  in  (3.3)  t  <  \z\ .  Then  the  rectangle  A  (z,  t)  is  located  below  {r  =  |£|}  and  above 
{t  =  0}  .  Hence,  we  obtain  from  (3.3) 


t—T~\~Z 


=  - 


dr  /  (A h,x,yui:j  +  dijUij)  (£,  r)  for  t  <  \z\  . 


T  —  t+Z 


Iterating,  we  obtain  from  here,  similarly  with  (3.4)  that 

(Ct)n 

\uij(z,t)\  <  — —  ,n=  1,2,... 

Hence,  uhJ  ( z,t )  =  0  for  t  <  \z\  .  The  same  way  uniqueness  of  the  problem  (3.3)  can  be 
proven.  □ 

Because  of  this  theorem,  we  can  consider  functions  Uij  (z,  t )  only  above  the  characteristic 
line  {t  —  |z|}  in  the  (z,  t)  plane.  Hence,  consider  new  functions  Wij  (z,  t )  =  Uij  (z,t  +  z) ,  z  > 
0.  The  domain  {t  >  z,z  >  0}  becomes  now  {t  >  0,z  >  0} .  Using  (2.4)  and  (2.6)-(2.8),  we 
obtain 

d2zwid  -  2dzdtwitj  =  -A h^yWij  +  ahJ  (z)  witj,  (z,  t)  e  {t  >  0,  z  >  0}  ,  (3.5) 

«4,j(u0)  =  ^,  (3.6) 

Wij  (0,  t)  =  Tij  (t) ,  dzWij  (0,  t)  =  Sij  (t) ,  t  G  (0,  T) ,  (3.7) 

Wi  j  G  C3  (z,  t  >  0) .  (3.8) 


4  The  Carleman  Estimate 

Consider  parameters  a  G  (0, 1/2)  ,/3,  v  >  0  as  well  as  a  sufficiently  large  parameter  A  >  1. 
We  will  choose  A  later.  Consider  the  functions  ^  (z,  t) ,  tp  (z,  t ) , 

ip  (z,t)  =  z  +  at  +  1,  ip  (z,  t)  =  exp  (A 'ip~u)  .  (4.1) 
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Define  the  domain  Dp  as 

Dp  =  {(z,t)  :  z,t  >  0,-0  (z,  t)  <  1  +  0}  .  (4.2) 

Then 

dDp  =  U tADp,  (4.3) 

d1Dp  =  {t  =  O,ze(O,0)},  (4.4) 

d2Dp  =  jz  =  0,  0  <  t  <  ,  (4.5) 

d3Dp  =  {z,t  >  0,if(z,  t)  =  1  +  P},  (4.6) 

\d3Df,=  exp  (A  (1  +  P)~v)  =min  <p(z,t).  (4.7) 

D/3 

When  applying  the  Carleman  estimate  of  Lemma  1,  we  will  have  Dirichlet  and  Neumann 
data  at  diD$,  as  it  follows  from  (3.7).  At  d3Dp  the  function  p  (z,  t )  attains  its  minimal  value, 
which  is  one  of  key  points  of  any  Carleman  estimate.  As  to  d\Dp,  we  will  not  have  any  data 
at  d\  Dp  when  applying  Lemma  1.  Also,  d\ Dp  is  not  a  level  curve  of  the  function  <p(z,t) . 
Therefore,  to  make  our  Carleman  estimate  valuable,  we  should  prove  that  the  integral  over 
d\ Dp,  which  occurs  due  to  the  Gauss  formula,  is  non-negative,  see  the  second  line  of  (4.8). 
The  latter  is  the  main  new  feature  of  Lemma  1. 

Lemma  1  (Carleman  estimate).  Let  a  G  (0, 1/2)  and  f3,v  >  0.  Then  there  exist  constants 
Ao  =  Ao  (ct,  j3,u)  >  1  ,C  —  C  (a,  /3,u)  >  0  such  that  the  following  Carleman  estimate  holds 


(uzz  —  2 uzt)2  p2dzdt  >  C A  /  {u2  +  u2  +  A 2u2)  p2dzdt 


Da 


Da 


+CA  J  [u2  +  \2u2)  (z,  0)  p2  (z,  0)  dz  (4.8) 

di  Dp 

— CA3  exp  [2A  (P  +  l)_l/]  J  (u2  +  u2  +  u 2)  dS, 

OsDf, 

VA  >  Ao,Vu  G  C2  ( Dp )  fl  {u  :  u  \q2 dp=  dzu  \d2Df,—  0}  .  (4.9) 

Proof.  In  this  proof  C  =  C  (a,  ft,  u)  >  0  denotes  different  positive  constants.  Consider 
a  new  function  v  =  up  and  express  uzz  —  2 uzt  via  v.  We  have 

u  =  v  exp  (—A , 
uz  =  (yz  +  \vif-v~1v)  exp  (— , 

vzz  +  2Xuif~v~1vz  +  X2u2if~2v~2  (l  —  — ~'tpv\  v  exp  [-Xif~v)  , 


U  zz 


Uzt  = 


Xu 


vzt  +  aXuif  v  1vz  +  Xuif  v  Lvt  +  aX2uz,if 


-1/-1 


1 2.  2  / — 2i/— 2 


exp  (-A^  v)  , 
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Uzz  -  2 uzt 


=  vzz  -  2vzt  +  (1  -  2a)  A  2v2ip 


2,  ,2./,— 2f— 2 


V+l) 


exp  (— A?/;-1') 


Hence, 


+  (2  (1  —  a)  Xvip  "  V2  —  2Az/?/>  "  Vt)  exp  (— . 


(uzz  -  2 uztf  tp2  >  2y2yi  -  2y3yx, 


(4.10) 


Vi  =  Vzz  -  2vzt  +  (1  -  2a)  A  v  ip 


2  2/ -2v-2 


y2  =  2  (1  -  a)  A vip  v  1vz,  y3  =  2Xuip  v  lvt. 


We  have 


2y2y\  =4(1  -  a)  Xvip  v  1vz  vzz  -  2vzt  +  (1  -  2a)  X 2v~ip 


i  22/ ,—2^—2 


1-^)» 


=  <92  (2  (1  -  a)  Xvip-^vl)  +  2(1 -  a)  Xu  (v  +  1)  ip~v~2v2 
+dt  (—4  (1  —  a)  Xvip~v~lv2P)  —  4a  (1  —  a)  Xv  (v  +  1)  ip~u~ 

i  a  \3,.3„/.-3i/-3  ft  V  1) „/,iA  „,2 


— 2,,2 


+d2  2  (1  —  a)  (1  —  2a)  X6v6ip 


+6  (1  —  a)  (1  —  2a)  A'V3  (z/  +  1)  ip~ 


1  -  v  v 

Xu  ) 


Thus, 


22/22/1  =  2  (1  -  a)  (1  -  2a)  Xu(u  +  l)ip  v  2v2z 

+6  (1  -  a)  (1  -  2a)  AV3  (u  +  1)  iP~3u~A  ( 1  -  (2t/  +  3)^A  v 2 

V  3Az/  / 


(4.11) 


+V  (-4  (1  -  a)  Xuip~u~1vf) 


+dz  2  (1  -  a)  Az 'ip~v~lv2z  +  2  (1  -  a)  (1  -  2a)  AV3^ 


^  U2  . 


Next  we  estimate  —2y3yi, 


21/32/1  =  -^Xuip  "  Vt  v22  -  2v2t  +  (1  -  2a)  AV2^ 


2,  ,2„/,— 2i/— 2 


("+!) 


1  - h±lV)  „ 

Az/  '  / 


=  <92  (— AXuip  "  Vtv2)  +  AXisip  "  V2tv2  —  4Az/  (z/  +  1)  ip  v  2vtvz 
+dz  (4Az/^-"-1n2)  +  4Az/  (z/  +  1)  ^"""V2 


3,,3„/,-3i/-3 


+<9*  —2  (1  —  2a)  X6vAip 


-6a  (1  —  2a)  A'V3  (z/  +  1)  ip~ 


1  -  h±li 

Az/  '  / 


Next, 

Hence, 


4A u^~u~1vztvz  =  dt  (2A uil)-v-1vfj  +  2aXu  (u  +  1)  ^~v~2v\. 
— 22/32/1  =  2Xv  (//  +  1)  (oiv2z  -  2vtvz  +  2v 2) 


-6a  (1  —  2a)  A3z/3  (y  +  1)  ip 


i  -3i/-4 


(2z/  +  3)  i  2 


3Az/ 


-ij)v  v 


(4.12) 


+dt 


2Xu^-v-1v2z  -  2  (1  -  2a)  A3z/V 


3,  3f— 3 


1  _  ++0  ) 
Az/ 


+.9,  [-4Az/'z/+I/_1'i+uz  +  4Az/'z/’_l'_1v2]  . 
Summing  up  (4.11)  and  (4.12)  and  taking  into  account  (4.11),  we  obtain 

(' uzz  -  2 uztf  +2  >  2y2yi  -  2 y3yt 

=  2Xu  (y  +  1)  ip~u~2  [(l  —  2a  +  3a2)  v2  —  2vtvz  +  2v 2] 

i  -  <*+V)  v> 


+6  (1  -  2a)2  AV  (z/  +  1)  ij} 


l-Zv-A 


3Az/ 


(4.13) 


+dt 


-2  (1  -  2a)  Az+/rI/"1u2  -  2  (1  -  2a)  A3z/ty 


3,  ,3„i.—3u—3 


i  -  k±lV )  t-2 

Xu 


+dz 


2  (1  —  a)  A uip  v  v2  +  2  (1  —  a)  (1  —  2a)  A3z/3,0 


3,,3„/ -3i/-3 


i  -  k±!V )  «2 

Az/ 


+<92  [-4Az/-0  "  Vu*  +  4Az/-0  v  1v2\  . 

Obviously  there  exists  a  constant  C\  =  C\  (a)  >  0  such  that 

(l  —  2a  +  3a2)  a2  —  2 ab  +  2 b2  >  C\  (a2  +  b2)  ,  Va  G  ^0,  ^  ,  Va,  b  G 

Hence,  integrating  (4.13)  over  Dp,  we  obtain 

[  (uzz  —  2 uzt)2  (p2dzdt  >  2Az/  (z/  +  1)  Ci  f  (u2  +  u2)  ip~u~2dzdt 


Da 


Da 


+6  (1  —  2a)2  A 3z/3  (z/  +  1)  /  i/j 


— 3f— 4 


1  - 


(2z/  +  3)  >  2 


3Az/ 


if;1'  v^dzdt 


Da 


di  Dp 


2(1-  2a)  Ai+T  u2  +  2(1-  2a)  A3z/ty 


3,,3„/,-3i/-3 


(4-14) 

Xu 
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cos  (n,  t )  dS 


+  / 

d3D, 

+  / 

dsDp 

+  J  [— 4A v/il)~l'~1vtvz  +  cos  (n,  z)  dS. 

dsDp 

Here  cos  (n,  t )  and  cos  (n,  z )  are  cosines  of  angles  between  the  unit  outward  normal  vector  n 
at  djDg  and  positive  directions  of  t  and  z  axis  respectively.  Since  the  number  v  >  0  is  fixed, 
we  can  incorporate  it  in  the  constant  C.  Change  variables  back  in  (4.14)  replacing  v  with 
u  =  vip.  Then  for  sufficiently  large  A  >  Ao  {v,/3)  (4.14)  implies  (4.8).  □ 

5  Proof  of  Theorem  2 

It  is  sufficient  to  consider  the  case  z  >  0,  since  the  case  z  <  0  is  similar.  Assume 
that  there  exist  two  pairs  of  vector  functions  ( u1,h(z,t )  ,a1,h(z ))  ,  ( u2,h  (z,  t ) ,  a2,h  (z))  sat¬ 
isfying  conditions  of  this  theorem.  Then  for  z,t  >  0  there  exist  two  pairs  of  functions 
( w1,h  (z,  t ) ,  a1,h  (^))  ,  ( w2,h  (z,  t ) ,  a2,h  (z))  ,  where  w1,h  (z,  t )  =  u1,h  (z,  t  +  z) ,  w2,h  (z,  t)  =  u2,h  (z,  t  +  z) . 
Denote 

wh  (z,  t)  =  w 1)h  (z,  t )  -  W2’h  (z,  t)  =  {wid  (z,  t)} (JJ)=(0i0)  * 
dh  (z)  =  a1,h  (z)  -  a2'h  (z)  =  {dij  (^)}jjj)f(0,0) 

Then  (3.5)-(3.7)  imply  that 

whzz  —  2 wzt  =  —Ah'Xtyivh  +  a1,h  {z)  wh  +  dh  (z)  w2,h  (z,  t ) ,  (z,  t)  G  {t  >  0,  z  >  0}  ,  (5.1) 


wh(z,  0)  =  0,  (5.2) 

wh  (0,  t)  =  0,  dzwh  (0,  t)  =  0,  t  G  (0,  T ) .  (5.3) 

By  (3.8)  wh  G  C3  (M  x  [0,  T]) .  Hence,  setting  in  (5.1)  t  —  0  and  using  (5.2),  we  obtain 

dh  (z)  =  -Adzdtwh  (z,  0) .  (5.4) 

Let  vh  (z,  t )  =  dtwh  (z,  t ) ,  v2,h  (z,  t )  =  dtw2,h  (z,  t) .  Differentiating  (5.1)  with  respect  to  t  and 
using  (5.3)  and  (5.4),  we  obtain  for  (z,t)  G  {t  >  0,  z  >  0} 

dhzz  -  2vzt  =  -  A htXtyvh  +  a1,h  (z)  dh  -  Adzvh  {z,  0)  v2’h  {z,  t ) ,  (5.6) 

vh  (0,  t)  =  0,  dzvh  (0,  t)  =  0,  t  G  (0,  T) .  (5.7) 
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-2  (1  -  2a)  \vi\)-v-xv2z  -  2  (1  -  2a)  A3i/V_3i'-3 
2  (1  -  a)  \vi)-v~xv2z  +  2  (1  -  a)  (1  -  2a)  A3zAr3l/-3 


+  1)  ,v  \  2 

1  —  - — 4 — )  v2 1  cos  (nt  z )  dS 
Xu 


Since  T  >  2 R,  then  (0 ,R/T)  C  (0,1/2).  In  (4.1)  choose  an  arbitrary  a  G  (0 ,R/T). 
Next,  set  in  (4.2)  /3  R  and  choose  an  arbitrary  v  >  0  in  (4.1).  Consider  equation  (5.6) 
for  vhi  rj  for  an  arbitrary  pair  (i,j)  G  [0,  Ni]  x  [0,  iV2] .  Next,  square  both  sides  of  the  latter 
equation,  multiply  by  the  function  p2  (z,  t )  and  integrate  over  DR.  We  obtain  with  a  constant 

M  —  M  (ji0,  ||a1,ft||C[0  ,  II^’^IIcOdr))  >  0  depending  on  listed  parameters 

J  (d^Vij  —  2 dzdtVij)2  (p2dzdt  (5.8) 

Dr 


<M  \yh  (z,  t )  2  p2dzdt  +  M  [v^j  (z,  0)]2  p2dzdt. 


D , 


D , 


Since  the  function  p2  ( z,t )  is  decreasing  with  respect  to  t,  we  obtain  from  (5.8)  and  (4.5) 


J  (d%Vij  -  2 dzdtVij)2  p2dzdt 

Dr 


(5.9) 


<M  \yh  (z,t)  p2dzdt  +  Mi  /  [vij  (z,  0)]2  p2  (z,  0)  dz, 


D, 


9iDr 


where  the  constant  Mx  —  Mi  (. M,R ) . 

Applying  Lemma  1  to  the  left  hand  side  of  (5.9)  and  using  (5.7),  we  obtain 


CA 


/  [(W  +  (W  +  A2PWV** 

Dr 


+CA  j  [( dzVij)2  +  A2  (Vij)2]  (z,  0)  p2  (z,  0)  dz 

9iDr 

-C  A3  exp  (2A  (R  +  l)-*')  J  [(9 zvitj)2  +  (dtvid)2  +  {v^f]  dS 

83  Dp 

<M  J  [vh  (z,t)]2  p2dzdt  +  Mi  J  [vij  (z,  0)]2  p2  (z,  0)  dz. 

Dr  9 1  Dr 

Choose  a  sufficiently  large  number  A0  >  1  such  that 

CA3 

max  (M,  Mi)  < 

Then  we  obtain  from  the  latter  estimate 

CX  j  [(dzVij)2  +  (d tVij)2  +  A2  (vij)2]  p2 dzdt 

Dr 


(5.10) 
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(Vij)2(z,  0)p2  (z,  0)  dz 


+CX3  J 

di  Dr 

- cx 3  exp  [2 A  (R  +  l)-*']  J  [( dzvi:j  f  +  {d tvitj)2  +  (vij)2]  dS 

d3Dp 

<M  J  \vh  (z,  £)] 2  (p2dzdt. 

Dr 

Summing  up  these  estimates  with  respect  to  (i,j)  G  [0,Ai]  x  [0, -A2]  and  using  (5.10),  we 
obtain  a  stronger  estimate 

J  (vh)2  (f2 dzdt  <  Cexp  [2 A  (R  +  l)-"]  J  (vz)2  +  (uf)2  +  (vh)2  dS.  (5.11) 

Dr  dzDp 

Let  £  G  (0,  R)  be  an  arbitrary  number.  By  (4.1)  and  (4.2)  p2  (z,  t )  >  exp  [2A  (R  +  1  —  e)~v] 
in  D r_£  and  Dr_£  C  Dr.  Hence,  making  the  estimate  (5.11)  stronger,  we  obtain 

exp  [2A  (R  +  1  —  s)~u]  J  (vh)2  dzdt  <  Cexp  [2 A  (R  +  l)~v]  J  (u^)2  +  (u^1)2  +  (u^)2  dS. 

Dr-6  dzDp 

Or 

J  (yh)2  dzdt  <  C exp  {— 2A  [(R  +  1  —  £)~P  —  (R+  l)-^]  }  J  {vz) 2  +  (u? ) 2  +  (vh) 2  dS. 

Dr-e  d'3  Dp 

Setting  here  A  — *  oo,  we  obtain 

f  (yh)2  dzdt  =  0.  (5-12) 

Dr-e 

Since  £  G  (0,  R)  is  an  arbitrary  number,  then  (5.12)  implies  that  vh  ( z,t )  =  0  in  Dr.  Since 
by  (5.4)  dh  (z)  =  —Adzdtwh  (z,  0)  =  —4dzvh  (z,  0) ,  then  dh  (. z )  =  0  for  z  G  (0,  R ) .  Finally 
since  by  (5.2) 


t 


o 


then  wh  (z,  t)  =  0  in  Dr.  □ 
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